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Abstract There is '^P2 neutron superfluid region in NS (neutron star) interior For a rotat- 
ing NS, the ■^P2 superfluid region is like a system of rotating magnetic dipoles. It will give 
out electromagnetic radiation, which may provides a new heating mechanism of NSs. This 
heating mechanism plus some cooling agent may give sound explanation to NS glitches. 
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1 INTRODUCTION 

Since the discovery of NSs (neutron star) and glitches in late 1960s, it is generally believed that there is 
superfluidity in NS interior (Ruderman 1976, Shapiro et al. 1983, Elgar0y et al. 1996). Cooling mechanism 
associated with superfluidity was first proposed by Flowers et al. (1976). Not until recently, the importance 
of superfluidity is considered seriously in the "Minimal model" (Gusakov et al. 2004, Page et al. 2004, 
Kaminker et al. 2006). 

Then one may ask: since cooling agent associated with superfluidity must be considered in the "Minimal 
model", what about the heating mechanisms? Heating mechanism accompanied with superfluidity has been 
discussed by Alpar et al. (1989), and was taken into NS cooling model by Umeda et al (1994), Page et 
al. (2006), Tsuruta (2006). Here we investigate another possibility: microscope magnetic dipole radiation 
(MMDR) in NSs, possibly another heating mechanism associated with superfluidity. 

The idea is as follows: there is ^P-z superfluid region in the interior of a NS. We calculate its paramag- 
netic properties in the presence of a background magnetic field. Since the neutron superfluid is in a vortex 
state, the ■^P2 superfluid region is like a system of rotating magnetic dipoles. This system will give magnetic 
dipole radiation. The emitted photon can not penetrate the NS matter, thus provides a heating mechanism of 
NSs. The origin of this heating mechanism is microscopic magnetic dipole radiation, so we call it MMDR 
heating of NSs. 

This heating mechanism plus some cooling agent may give sound explanation to NS glitches (Bai et al. 
2006; Peng et al. 2006; Peng 2007). 

Before going into detailed calculations of MMDR, we will first look at the superfluidity in NSs. 

2 SUPERFLUIDITY IN NS 

There are two relevant regimes of neutron superfluid inside the NS interior: one is the isotropic ^5*0 
neutron superfluid within a density range 1 x 10^'^ < P {g cm^^) < 1.6 x lO^"*. The critical temperature is 

^,(l5o)«lxlolOA^ 

Another important regime is the anisotropic '^Pj neutron superfluid within a wide density range 1.3 x 
10^* < p (gcm^^) < 7.2 x 10^"'. The critical temperature is 

T,(3P2) = A,„,,(3p2)/2fc « 2.78 X 10*^ K. (1) 
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We note that the energy gap A{^P2) is almost a constant about the maximum with an error less than 3% in 
a rather wide density range 3.3 x 10^^ < P {g cm~^) < 5.2 x 10^^ (see Fig. 2 of Elgar0y et al. 1996, but 
we neglect the F state of neutron Cooper pair here). 

It is well known that a rotational superfluid must be in the superfluid vortexes. In general the vortex 
filaments are arranged in a symmetric lattice, they are parallel to the axis of rotation of NS almost rigidly. 
The circulation of every vortex filament F is quantized 

F = jv-dl^ nFo, (2) 

Fo = (3) 
2m„ " 

Here n is a circulation quantum number of the vortex, m„ is neutron mass, h is the Planck's constant, Fq is 
the unit vortex quanta(Feynman 1955, Lifshitz et al. 1999). 

It is supposed that the core of the superfluid vortex is a cylindrical region of normal neutron fluid 
immersed in the superfluid neutron sea. As usual, the radius of the core of vortex ao, is taken to be the 
coherent length of neutron superfluid 

Here Ep is the Fermi energy of the neutrons, fc^? is the corresponding Fermi wave number, p is the total 
density of the neutron superfluid region(Ruderman 1976). Outside the core of the vortex, neutrons are in a 
superfluid state. 

The kinetic properties of superfluid vortexes follow directly from the Feynman circulation theorem 
eqlU The superfluid neutrons revolve around the vortex line with a velocity 

Vs{r) = , (5) 

where r is the distance from the field point to the axis of the vortex filament. The distribution of the angular 
velocity of the neutron superfluid revolving around the vortex filament is 

^s{r) = (6) 

Therefore the revolution of the superfluid neutrons around the vortex filament is placed in a differential 
state: the closer to the center of the vortex, the faster the superfluid neutrons rotate. Near r « ap, the 
angular velocity reaches its maximum value 

^ ^ nh 
2TO„ag 

And inside the core of the vortex r < ao, the normal neutron fluid revolve rigidly at angular velocity of lUc- 
According to Feynman (1955), the number of superfluid vortex filament per unit area is 2ri/nFo, where 
51 is macroscopic angular velocity, n is mean circulation quantum number Then the average separation h 
between vortex filaments and the total number of the superfluid vortexes in the interior of NS are respec- 
tively 



1/2 



iV„..= (- =-^i?^. (9) 



Here R is the radius of the superfluid region. 
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In this paper, we only consider the thermodynamic condition. That is n = n = 1. Although deviation 
from equilibrium have been discussed by Yakovleve et al. (2001), Reisenegger et al.(1995, 2006). 

For a neutron superfluid vortex in the interior of a NS, ao ^ 10~^° cm and6^ 10~^ cm. Therefore 
the vortex core is very tiny, the distribution of the vortex filament is exceedingly sparse. The separation 
between vortex filaments reaches a macroscopic scale, and superfluidity is just a macroscopic quantum 
phenomenon. 

The magnetic behavior of superfluid region will be discussed in the next section. 



3 INDUCED PARAMAGNETIC MOMENT OF THE NEUTRON SUPERFLUID IN THE 
B-PHASE 

3.1 Two phases of the neutron superfluid 

A neutron Cooper pair has spin angular momentum with a spin quantum number, s = 1. The 
magnetic moment of the '''P2 neutron Cooper pair is twice that of a neutron, 2/x„ in magnitude, where 
^,n = —0.966 X 10"^'^ erg gauss~^ is the anomalous neutron magnetic moment. And its projection on an 
external magnetic field (z-direction) is x 2 /i„, = 1, 0, —1. It is interesting to note that the behavior of 
the ^P2 neutron superfluid is very similar to that of the liquid ^He at very low temperature (Leggett 1975): 

1. The projection distribution for the magnetic moment of the ^P2 neutron Cooper pairs in the absence of 
external magnetic field is stochastic, or "Equal Spin Pair" (ESP) phase. The ^P2 neutron superfluid is 
basically isotropic and with no significant magnetic moment in the absence of external magnetic field. 
We name it as the A-phase of the ■^P2 neutron superfluid similar to the A-phase of the liquid ^He at 
very low temperature (Leggett 1975). 

2. However, the projection distribution for the magnetic moment of the ^P2 neutron Cooper pairs in the 
presence of external magnetic fleld is not stochastic. The number of the ^P2 neutron Cooper pair with 
paramagnetic moment is more than the ones with the diamagnetic moment. Therefore, the ^P2 neutron 
superfluid has a net induced paramagnetic moment and its behavior is anisotropic in the presence of 
external magnetic fleld. We name it as the B-phase of the ^P2 neutron superfluid similar to the B-phase 
of the liquid ^He at very low temperature (Leggett 1975). 



3.2 Induced paramagnetic moment of the ^P2 neutron superfluid in the B-phase 

We now consider the paramagnetism of ^P2 neutron superfluid regions. FoUowing standard treatment 
of magnetism (Pathria 2003, Feng et al 2005), the Hamiltonian of the system in the presence of external 
fleld is 

H = -2/2„ • B = -2^ir^zB. (10) 

Here B is the external fleld (in the z direction), 2^„ is the magnetic moment of the ■^P2 neutron Cooper 
pair, 2/u„z is its projection on the z direction. 

The ensemble average only gives an additional thermal factor 

(2m„) =2Ai„/(MnB/A;T) (11) 

fil^nB k T) = ^ J": . (12) 
1 + 2 cosn/?2/x„i3 



f{^^B/kT) = ^^. (13) 



In the limit of high temperature 

The oc y behavior is just the Curie's law of paramagnetism in terrestrial laboratory. The qualitative behavior 
is : as a NS cools down, its internal temperature T decrease, while the thermal factor increase. The '^P2 
neutron Cooper tends to align in the same direction. This is the mathematica formalism of the B-phase of 
the ^P2 neutron superfluid. 
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As mentioned by Lifshitz et al. (1999), there is a finite probability for two neutrons to combine into a 
Cooper pair Since only particles in the vicinity of the Fermi surface contributes, only a finite fraction q of 
the Fermi sphere are in the condensate state 




(14) 



~ 0.087. 



We have used the relation: Afc = ^/2mnA,pF = \/2mnEF- Here A ^ 0.05 MeV is the energy gap of the 
corresponding superfluid region, Ep ~ 60 MeV is the Fermi energy of the neutron system, pi? is the 
corresponding Fermi momentum, AA: is the thickness of the shell which will combine to Cooper pairs. 

In conclusion, for a specific volume Ay of •^P2 neutron superfluid region, only a small fraction of the 
Fermi sphere can combine into Cooper pairs. Ensemble average gives another thermal factor. So the 
net magnetic moment of the specific volume is 

Mav = AVpNAfin qfitiuB/kT). (15) 

Here Na is Avogadro's constant. The volume /S.V is specially chosen so that it is large enough to contain 
an numerous number of Cooper pairs, but small enough compared with macroscopic scale. It is a meso- 
scopic volume. This is required by the coherence calculation of the microscopic magnetic dipole radiation 
(MMDR) below. 

For a rotating NS, the superfluid region is in vortex state. It have paramagnetic moment in the 
presence of external field, on the mesoscopic scale. So the neutron superfluid vortex is like a system of 
rotating magnetic dipoles, which will give magnetic dipole radiation. This radiation can not penetrate the 
NS matter, thus provides a heating mechanism of NS associate with superfluidity. 

The detailed calculation of MMDR is the scope of the next section. 

4 MICROSCOPIC MAGNETIC DIPOLE RADIATION HEATING 
4.1 A working assumption 

In a superpluid vortex, each superfluid neutron revolve round the axes of the vortex with angular ve- 
locity Lo{r) (eq|6]l. It is well known that a rotating magnetic dipole will give out magnetic dipole radiation 
(MDR). Therefore, the superfluid neutron will give electromagnetic radiation. The frequency of emitted 
photons is equal to the rotational frequency of the neutron. 

We will use a phenomenological method to explore the radiation problem in this paper (see, e.g. 
Feynman 1955; Androdikashvili et al. 1966). The process is as follows: the rotational velocities of su- 
perfluid neutron will decrease during the emission of the magnetic dipole radiation due to the dissipation of 
its energy. The superfluid neutrons, hence, will drift out according to eq|6l and then the transverse pressure 
exerting on the normal neutrons in the vortex cores will decrease. The normal neutrons will move out to 
r > ao and become superfluid ones. At the end of the vortexes, other normal neutrons located in the normal 
neutron layer in the interior of neutron stars will flow into vortex cores along the axes. At the same time 
the densities of superfluid neutrons at the boundaries of vortex lattices will increase driving a flow into the 
normal neutron layer along the direction perpendicular to the axes. A "local circulation" will be formed in 
the superfluid vortex region in this way. This process is very similar to the Ekman pumping (Greenspan 
1968, Anderson et al. 1978). The normal neutrons in the crust of the neutron star will drift inward during 
the process of these local circulations, and the crust will slowly shrink a little. We may suppose that the 
energy of the magnetic dipole radiation is really transformed by the released gravitational energy of the 
inward moved crust through the local circulation. 
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4.2 Microscopic magnetic dipole radiation heating 

First, we calculate the magnetic dipole radiation by one single vortex. The power radiated by one 
neutron is (Huang et al. 1982) 

/ 

o, .4 

= ^T.^\Mi\f){f\mf) 

2^4 

= —(^\Ml^'^^). (17) 

where c is the speed of light, w is the angular velocity of superfluid neutron, |/) is the final state, \i) is 
the initial state, and is the operator of magnetic moment. If we consider a coherent small volume At/, 
and take 

{i\MlM^\i) = {i\Ml\i) « Mly sin^ 0. (18) 

Here Mav is the corresponding paramagnetic moment in a specific volume of Ay, 6 is the angle between 
the background field B and the rotational axis VI. Thus the required magnetic dipole radiation rate is 

VFav = ^IMavP sin^e, (19) 

The formula is similar to the classical case (Shapiro et al. 1983), because we have made approximations in 
obtaining egfTSl 

In calculation of the radiation power, the coherent effect must be taken into consideration (Peng et al. 
1982). This is like the coherent effect in electromagnetism. We first calculate the contribution of a specific 
volume ISV, then make a summation of the all the specific volumes. This is done for one superfluid vortex. 
The total power of MMDR is simply the sum of all single vortex contributions. See the appendix for details. 

The total power of microscopic magnetic dipole radiation (MMDR) is 

Wtot = Nyort ry log — . (20) 

Here Nyort is total number of superfluid vortexes eq|9] rj is the efficiency of coherence. Radiation of one 
single vortex is loe; — , where the contribution factor A is 

o 4 

sin2 e\pNAPinqf?Rojlal. (21) 

The factor A is proportional to R, the vortex number N^ort is proportional to R^, as can be seen in eq|9] So 
the total power Wtot is proportional to R^, that is proportional to the volume of the neutron superfluid 
region. This is what it should be. 

The essential and discussion of MMDR heating is presented in the next two sections. 

5 MMDR VS. OTHER HEATING MECHANISMS 

Heating mechanisms in NSs can be classified into different categories according to their energy input. 
There are several kinds of energy input: magnetic, rotational, chemical, and confinement energy etc. The 
corresponding heating mechanisms are respectively: 

1. Ohm heating (Page et al. 2006 and reference therein); 

2. Vortex creep heating (Alpar et al. 1989, Umeda et al. 1994); 

3. Retochemical heating (Reisenegger 1995, Reisenegger et al. 2006). 

4. Deconfinement heating (Yuan et al. 1999, Kang et al. 2007); 

MMDR heating has two distinguished points compared with the heating mechanisms stated above: 

1 . It is a heating mechanism associated with superfluidity, thus no superfluid supression. 

2. Its energy input is gravitational energy as stated in our working assumption. 
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6 CONCLUSIONS 

Here we have presented another possible heating mechanism of NSs, associated with superfluidity. It 
can be compared with other heating mechanisms and cooling agents of NSs (Gusakov et al. 2004, Page et 
al. 2006). 

HMDR heating VS cooling 
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Fig. 1 MMDR Heating VS Cooling. The brown one represents MMDR heating. It increase with 
decreasing temperature, as a result of increasing thermal factor The red, green, and black one 
corresponds to photon cooling, MUrca (Modified Urea) process, and PBF (Pair Breaking and 
Formation) process respectively. (Adapted from Gusakov et al. (2004).) 

As shown in Fig. [T] MMDR heating may be dominate only in the photon cooling stage. So it will not 
affect the cooling scenario of young and mid age NSs. But for old NSs, e.g. PSR 1055-52, it may serve as 
a moderate heating (Tsurata 2006, Page 2006). 

Using toy model given by Yakovlev et al. (2003), we can make illustrative calculation of NS cooling, 
including the MMDR heating. Fig. |2] shows that, there is a plateau in the late stage of NS cooling curves. 
This plateau can be compared with that of Kang's (2007). 

This MMDR heating must cease after some time. In our case, the cease of MMDR heating has several 
possibilities. 

1. As shown by Huang et al. (1982), the energy input of MMDR heating is through Ekman pumping. 
When this agent is out of work, MMDR heating have to make a pause. 

2. When there is a phase transition in the core, e.g. deconfinement of hardrons, all the hadron processes 
disappear including MMDR heating. 

3. Due complexities of superfluidity gap (Elagr0y et al. 1996), the •^P2 superfluidity region becomes slimer 
when the core becomes more compact. The MMDR heating contribution is ignorable in this case. 

We presented another possible NS heating mechanism associated with superfluidity. The exact effect 
of MMDR heating needs accurate calculation of the cooling curves. This heating mechanism plus some 
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Fig. 2 Cooling Curves Including MMDR Heating. The lower and upper curve correspond to 
MUrca and bremsstrahlung process dominated case respectively. Whereas the plateau at latter 
stage is due to MMDR heating. The physical input are the same as Yakovlev et al. (2003) except 
that we include the MMDR heating. Remind: these cooling curves are for illustrative use only. 



cooling agent may give sound explanation to NS glitches. The detailed investigation is the scope of another 
work (Peng 2007). 
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Appendix A: DETAILED CALCULATION OF MMDR POWER 

We will follow the routine of Peng et al. (1982). At radius r from the center of neutron superfluid 
vortex, the neutron rotational frequency is 

As stated in the main text, the frequency of radiated photons is equal to the neutron rotational frequency 
Ws{r). The wave length of the radiated photon is 

A. r = = (r > ao). (A.2) 

The length of one single superfluid vortex is H = ^R, where R is the radius of the superfluid region. One 
superfluid vortex can be separated into several segments, we need only consider the coherent effect inside 
each segment. If we assume the a length scale Ag (r), then a cylindrical shell between the radius r and r + dr 
can be cut into H/\g (r) segments. Here the specific volume AV is 

AV = 2'!rrdrXs{r). (A.3) 

Paramagnetic moment in the specific volume is 

Mav = AVpNAfinqf (^) ■ (A.4) 



This was prepared with the ChJAA macro vl.2. 
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Using eq[T9] differential power of one single superfluid vortex is 



dW^W^^-^lMAvpsin^e (A.5) 
= '^T^smH\pNAtinqf\^\27rrdrX,{r)\^ (A.6) 



3c3 

3c- 



H Xsir) sin' 9\pNAtinqfrr'dr' (A.7) 
^ sin' 9\pNA^inqf\''Ru^r^dr\ (A.8) 



In obtaining the final expression, we have used the definition of As(r) and H. Using non-dimensional 
quantities 

r = aor' ^^^^ 
, 1 

the differential power is 

dVF(i) = sin' e\pNA^iuqf?RLol4uJ'\"dr" (A. 10) 

= Au'^r"Ar" (A. 11) 

= AAI (A. 12) 

When integrating d/, a (5-function is included automatically, 5{lo' — 1/r"). Performing the integration from 
flQ to h gives i log The contribution factor A is 

o 4 

^ = — sin^ e\pNAPinqf\^Rojlal. (A.13) 

OC 

In the end, the microscopic magnetic dipole radiation power of one superfluid vortex is 

W^^^=A-\og—. (A. 14) 

2 ao 

In the above calculation, we assumed a length scale As(r). More generally, if the length scale is taken 
to be rjXs (r), where rj is an efficiency factor, the MMDR power of one superfluid vortex is 

W^^^ =ri A- log— . (A.15) 
2 ao 

Discussions in Sect|2]provides the total number of superfluid vortexes 

N,ort={^\ ^i?'. (A.16) 

In conclusion, the total power of MMDR in the neutron superfluid region in the interior of NS is 

Wtot = N^ortriA]- log — . (A. 17) 

2 ao 

This is what we use in the main text. 



